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■ A discussion of the overlap problem of reweighting approaches to evaluating critical 
CO ' phenomenon in fermionic systems is motivated by highlighting the divergence of 

' the joint probability density function of a general ratio. By identifying the bounds 

for which this integral can be expressed in closed form, we establish criteria for 

■ accurately mapping the joint ratio distribution of two disjoint ensembles through 
interpolation. The approach is applied to QCD with four staggered flavours to 

^ l' evaluate the critical line in the /3 — fi plane. 



X 



1. Introduction 



The pathology of fermionic reweighting schemes can be succinctly expressed 
5^ , in the free energy density difference, A/, given by the ratio of two partition 

functions for a given finite system of intensive and extensive variables, (3, 
and, V. Reweighting methods attempt to effect this difference through 
numerical simulation by normalising observables with the ratio of the Monte 
Carlo functional measures for the two separate ensembles. 

^ = exp{-Af PV} (1) 

However, as the free energy density difference that can be evaluated 
is necessarily positive the exponent of this ratio can become vanishingly 
small. The relevance of the ensemble that can be numerically evaluated, 
n', to the phase space of the ensemble of interest, il, therefore becomes 
questionable. Simple-mindedly we could ask how close two ensembles must 
be for a reweighting measurement to be unaffected by this finite difference. 
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2. Probability Density Function of a Ratio 

This question can be expressed in general terms through the joint prob- 
abihty density function, 4'A{t), of the ratio of two normally distributed 
variables, x, and, y. Both having a given mean, fj. , and variance, a, and 
in the following example below vanishing probability densities everywhere 
but for positive values of y and x. 

y 

t = - y ■■ 4'a{iJ'A,(^a) X : <j)B{lJ. 8,(^3) (2) 



X 



POO i'{t-\-at)x 

d(j)A{t)= / dx (j)A{tx).(j)B{x) / dy (3) 

J x—O •J y—tx 

d(j)A{t) 

\t) = lim — ^— = dxx (j)A{tx)4B{x) (4) 

' dt^o at 7^=0 

_ a\iiBt + cr%fj,A { {p-B - iJ'Atf 



Even for this overly simplistic case the joint probability density function 
of the ratio has a nonzero Cauchy component and so all higher moments 
(including the mean and variance) arc undefined. Insight is gained with the 
central limit theorem under Lyapunov conditions ^, where the distribution 
of the ratio will asymptotically approach normal if the mean of the ratio 
is several standard deviations from zero. Similarly, under these conditions 
the higher moments can thus be properly defined for the discrete finite free 
energy density distribution of Eq. (1) which relates to the pathology of 
reweighting ^. 

i i 

{50f_f{0,0} ^\ , ^{1,1} A 



(o)/ V {OV J \ {1}^ ; \{o}{i} 

Ignoring the integrated autocorrelation times of the Monte Carlo eval- 
uation in the above defintion the relative error of a reweighted observable 
(O)^ is clearly minimised for when the observable O approaches unity. 
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Conversely, if an observable is normalised so that the mean is unity with a 
small standard deviation, we have ensured that Lyapunov-type conditions 
are valid (at least for the numerator). Since with reweighting a measure- 
ment can be redefined relative to a different ensemble simply by explicitly 
evaluating a free energy density difference in the measurement, the remain- 
ing denominator can be expressed as the product of a series of terms each 
incrementally close to unity. It is then straightforward to show with Eq. 
(8) that the relative error of such a product converges with an increasing 
number of increments for the given O factored in this manner. By express- 
ing the free energy density difference of Eq. (1) as the product of series 
of vanishingly small differences under appropriate contraints the finite free 
energy density difference of reweighting can thus be interpolated. 



3. Sign Problem 

The importance sampling evaluation procedure of a Monte Carlo method 

is essentially undefined for non-positive definite weights. This is the case 
for spin systems such as the Hubbard model where one effective measure 
which is used for simulations is the modulus ensemble ^. 

_ jDUOdetMe-^ 

jDUdetMe-s ^''> 

detM is real, sgn{detM) = | (10) 

_ HdctMHe-^ 
^^^11 " JDU ||dctM||e-S ^^^> 
(O.S5n(dctM))|| 
" (s<7n(detM))|| ^^^^ 

To interpolate sgn{detM) we would therefore rewrite the ratio of expec- 
tations in Eq. (12) as a product of expectations of increments close to unity. 
Bringing the numerator into a similar product form as the denominator if 
required by a significant probability density of the numerator at zero. 

^ ' Oo - l=NS (13) 



Oo 



. s5n(detM) 



{O) = — ^ ^ (14) 

nti Re (i^^ .exp{^{l + sfln(detM)) } 
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4. Overlap Problem 

Following the inclusion of the chemical potential n into the fermionic ac- 
tion, lattice QCD is similarly unamenable to direct Monte Carlo treatment 
as detM(/i) is complex valued for ^ ^ 0. A case in point is the finite den- 
sity Glasgow method The distribution of a set of normalised expansion 
coefficients is wanted at a point on the critical line, /xi, but an ensemble 
can only be generated on the real line for /Xq = 0. We may now, though, ex- 
press the relation between the two normalised ensemble-averaged expansion 
terms for these two regions (/ii, /io) in terms analogous to the sign/modulus 
relation defined for the sign problem in Eq. (14). 

/ c„ dctM(pi) \ 
\detM(/Xl) /„ ~ / detM(Mi) \ ^ ' 

As before to accurately map the free energy density difference between 
ensembles a normalising factor is inserted to bring the numerator close to 
unity and the remaining free energy density difference is interpolated in the 
denominator. 



detM(^i) 



dctM(/io) detM(p2) 



Mo 



/ detM(/ti) 1 
\detM((Uo) detM(/H2) 



detM(/J3) 

/ dctM(/Ji) 
\detM(,*3)/^^ 



nN I dctM(nAjj) 
n=l \ detM( [n+l]A^l) 



(17) 



(18) 



detM(^i 



detM(/Xo)detM(/x2) ~ detM(/x3) ^ detMifis) " ® ^ ^ I J 

For convenience, rather than determine the normalised distributions at 
the critical line from ensembles generated for several /3 values we addition- 
ally reweight in The modified coefficient being again normalised to 
unity through a choice for /is. 

c„ ^ c„ . exp{-Sgi0i) + SgiM} (20) 

The zeros of this polynomial which approach the real axis of the ex- 
pansion variable in the thermodynamic limit identify the critical line ^. A 
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Figure 1. The critical line in the fSi — fx plane for ny = 4 QCD, determined from the 
zeros of a fugacity polynomial expansion. 

preliminary critical value for /i is plotted as a function of /3i in Figure 1. 
It should be noted that since the polynomial is deflated during rootfinding 
there is no actual dependence on the transition value of /^i , and consequen- 
tially no need to tunc the rcwcighting parameters to effect cancellations 
through the covariance. The above ensemble consists of 2,000 configura- 
tions on a 4^ volume at f3o = 5.04 with m = 0.10, though strictly the 
quadratic fit ^ is for a smaller bare mass (m = 0.05). The congruence is 
therefore qualitative, although the fall-off at /i ^ 0.7 ^ is perhaps more 
consistent with expectations for the density of nuclear matter at /3i = 0.0. 
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